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A generalized supersymmetric representation of the Hubbard operator al- 
gebra is considered. This representation is applied to the infinite-C/ Hub- 



^ ■ bard model. A mean-field theory which takes into account both on-site and 

^~J ■ inter-site virtual boson-fermion transitions is developed. Unlike previous ap- 

f~^ ■ proaches, the mean-field theory considered is free from divergences. A possible 

o. 

^D . application of these results to the ferromagnet-paramagnet transition, as well 



as to other problems is discussed. 



a 

O . I. INTRODUCTION 

'k>( I The description of magnetism within the Hubbard model is the long-standing problem, 

c^ I At f/ = oo the ferromagnetic (but not antiferromagnetic) state can occur, and, as shown by 

Nagaoka, for small hole concentrations 6 this is the ground state. The stability of saturated 
ferromagnetism for not too small 6 was investigated within different approaches (see, e.g. 
Refs. fTj-Q). The critical concentration of holes was estimated as 6c — 1/3 for the transition 
into non-saturated ferromagnetic state, and as 6'^ — 2/3 for that into non-magnetic one. The 
same result for 6c can be obtained by comparing the total energies in the slave-fermion ap- 
proach (which describes correctly the excitations near half-filling) and slave-boson approach 
(which is suitable for the non- magnetic phase), see Ref. P] and references therein. 

Slave-fermion and slave-boson approaches have the advantage that they correspond to the 



N -^ oo limit of the generalized U = oo Hubbard model with A^ fermion flavors (which will 
be referred to as SU{N\1) model), and 1/A^ corrections can be found in a regular way. These 
two approaches treat different kind of excitations, and therefore have different applicability 
regions. While slave-fermion approach is convenient at small hole concentrations and enables 
one to describe naturally the magnetic state, the slave-boson approach works well at high S. 

To obtain the physical picture at intermediate hole concentrations, a supersymmetric 
approach should be developed which interpolates between slave-fermion and slave-boson 
ones. The situation here is similar to that in s — / model where boson representation of 
impurity spin operators describes correctly magnetic phases, while fermion representation is 
suitable for description of Kondo (i.e. nonmagnetic) state. To describe continuous transition 
between these phases, the supersymmetric approach of Ref. can be used. 

In this paper we consider the application of the supersymmetric approach to the SU{N\ 1) 
model and obtain the effective action of the infinite-^ Hubbard model in supersymmetric 
representation. Then we consider the mean-field approach to this action and obtain the 
corresponding self-consistent equations. 

II. REPRESENTATIONS OF THE HUBBARD ALGEBRA 

The standard way of treating large-?/ Hubbard model is introducing the Hubbard oper- 
ators X°^ («,/5 = 0, ±). These satisfy at a lattice site the commutation relations 

[X"^ X^']± = dp^X"^' ± 5a5X^^ (1) 

and the constraint 

^X"" = l (2) 

a 

The operators X""^ give the spin S = 1/2 realization of the SU{2\1) superalgebra. To 
construct 1/N expansion we have to generalize this algebra to SU{N\1) by introducing the 
operators x"'^ {ci,P = 0...N) which satisfy the same commutation relations (|I|). We also use 
the generalized form of the constraint 



Ex"" = Qo 



(3) 



which reduces to standard one, (0), for the physical case Qq = 1. The operators x"" yield 
a representation of the SU{N\1) superalgebra with the "superspin" Qo/N. 

The slave-fermion representation of the operators x"^ through A^ bosons and one fermion 
has the form 



x'"'' = b%,,x 



o-O 



btf, /° = ff 



(4) 



{(J, a' = 1...N), with b„ and / being Bose and Fermi operators respectively. For Qq < N 
there exists also the slave-boson representation 



X"^' =cic^', X'^° = 4a, x'' = a^a 



(5) 



where c^ and a are Fermi and Bose operators. 

To obtain the supersymmetric representation which interpolates between slave-fermion 
and slave-boson ones, we introduce, following to Refs. [S-TO, the operators 



[a = 1...N) with 



^. 






\&r 



L\ 



\fj 



(6) 






(7) 



To make a distinction between the representations with different symmetry, we consider the 
second-order Casimir operator of SU{N\1) 



C2 = x"""" r "" - r^^"" + x^'^r'' - {x^y 



(8) 



Expressing this in terms of \1/ we obtain 



C2 = ^i^o\N - 1 - 4r3^^] - [6, et] 



(9) 



where T3 is the Pauh matrix, 



e = blc, - p 



(10) 



is the mixing fermion-boson operator with {0, 0^} = Qq, and the summation over repeated 
indices is assumed. Taking into account the constraint (0) we have finally 



C2 = Q,{N-l-Y) 



(11) 



where 



y = vi/tr3vi/^+ [e,et] 

By fixing Y = — Qo + 1 ••• Qo — 1, "we obtain representations with different symmetry. 



(12) 



III. SU{N\1) GENERALIZATION OF THE U = 00 HUBBARD MODEL 



The Hamiltonian of the SU{N\1) model can be now rewritten as 



^ = Et^JXfxT = Y.t^M.'^^0^]o^,a 



(13) 



For A^ = 2, Qo 
function reads 



1 this coincides with that of the U = 00 Hubbard model. The partition 



Z= [ D^expl- f dr{Co + H) 



(14) 



where 



2C 



-^0 = E ^L f I- + -^^ + ^'^-^ + ^^A "^^--itT. ©1©^ - E(0>^ + A.go + ^J^5) (15) 

is the free Lagrangian, 5 is the concentration of holes. Following to Ref. p], we perform the 
replacement ^ —* g^ with 
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1 + 7777/2 T] 

^ -V 1 - w/2 J 
4 



(16) 



to obtain the gauge-invariant Lagrangian in the form 



i 

vo i i 

Performing the Hubbard- Stratonovich transformation, we obtain 0] 



^o = E^ 



dr + Xi + flSaO + 



Ci D^a, ' 



V 



/ 



vl/. 

^ 7 



+Qo E c^lDoai - ElCii" + A^Qo + f^S) 



(17) 



(18) 



with Dq = dT- + 2Q. The Hamiltonian ([131) is decoupled in the same way as in Ref. |TT|, and 
we derive 

+'Pijaif] + /ia]Pji + c\j)jaQji + Qii&Lci<T 
= ^ U, \^lV,j^,. + *lo^i,-^,o - STr(V,,Z,,)l (19) 



»j 



where STr(...) is the supertrace; V^V, Q, Q are independent Grassmann variables, 






(20) 



and 



V. 



(aq^ 



7- 



V 



C — "P 
-V B 



U 



(21) 



y / 



The model (14) with (18) and ([T9|) is invariant under the gauge transformation 



a ^ a + (9^ + 2Ci)r7i 

l^jj *■ 9% '^ij9j ; ^ii ^ 9i-^ij9j 



(22) 



The action S = C'^ + 'H can be rewritten in the form 



where 



5 = E^ 



+E^I 



{dr + Xi + CiT3)6ij + _ 

^ Qij ^ijt J 

( 



^ 



jo- 



C-f P 



y -* «i ^ir J 



^r 



io 



+Qo E alDoai - STr(V,,Z,,) - J2{QY + \Qo + ^6) 



(23) 



Qji = Doai6ij + tQji, Q^^ = D^aidij + tQji 
Pji = Doai6ij - tVji, Pij = Doal6ij - tV 



«j 



(24) 



IV. MEAN-FIELD THEORY 



Our purpose now is to consider the mean- field approximation for the action (p3|), which 
does not violate the gauge invariance (p^). To this end, we will take into account the fiuctu- 
ations of Q,V and a, while other fields will be considered in the mean-field approximation. 
The motivation for this approximation is as follows. 

(i) The fields Q,V and a are equally important: while the field a describes the on-site 
virtual transitions of bosons into fermions and vice versa, the fields Q and V describe the 
same transitions with simultaneous intersite hopping. As can be seen from (p3|), Q and V 
play the role of "spatial components" of the gauge field, while a is only its time component. 

(ii) All the fields, besides Q, V and a, being taken into account at the mean-field level, are 
shifted properly by the gauge transformation. At the same time, the fields Q, V and a have 
zero mean-field value due to their fermionic nature and therefore can not be transformed at 
the mean- field level. 

Thus, Q, V and a make the minimal set of fields, fluctuations of which should be taken 
into account to keep gauge invariance. The gauge transformation of resulting theory is 



considered in Appendix. Note that taking into account only fluctuations of the field a leads 
to divergences due to violation of gauge invariance |T0[ . 

Unfortunately, unlike treating of single- impurity problem of Ref . |^ , the fields Q and V 
can not be completely removed by gauge transformation. Thus it will be more convenient 
for us to work in the gauge where a^ = 0. 

Integrating over \E^jq,, expanding the action to second order in Pjj, Qij at «« = we obtain 
the effective action in the form 

Seff = -iVSTrln[G,-i(q,z^„)G;i(q,2^„)] - STr ln[G-i(q,za;„)G'7i(q, za;„)] 



Cl,iu)n (5,(5' 






7ls'{q,iuJn) (25) 



i 

where TZs = (Q5,Vs), and the polarization operators are given by 

nf,,(q, toOn) = NJ2 GbiK zz/„)a(k + q, iPn + z^„)e^'^(^-^') 

k,iv„ 



NE' 



''^k + ""-k+q JkjS-S') 



-c I r-b 



e 



k 1'^n ^k+q + ^k 



k,JDn 



E 



^k-q + ""-k MS-S') 



J 



e 



k "^^n — £^k-q + ^k 

with G;,c,a,/(k, iz/n) are standard Bose (Fermi) Green functions of corresponding fields with 
the spectra 

£k = -Btk + A - C + /i 

^k = -^^k + A - C (26) 

The mean-field values JF, A, C and B have to be determined self-consistently 



C={clc,,), A={ala,) (27) 

For A, ( and fj, we have the constraint equations 

(///.) + {aU^) = 5 

{blha) + (cLq^) = 1-6 
(CU + (///.) + 7^(^M)= 25 (28) 

Further we consider the most interesting case of a "mixed" phase where both the bosons 
a, b are condensed. Physically, this corresponds to a non-saturated ferromagnetic state. The 
then concentrations 6c, 6'^ where the condensates of b and a bosons vanish will determine the 
transitions into saturated ferromagnetic and nonmagnetic states respectively. 

First, we consider the standard (noninteracting) mean-field theory which neglects the 
fluctuations of the Q and V fields, i.e. the second line of (pSf). We obtain in this case at 
T = 



na = Aq = 6 -Y^ 



n 



f 



n, = -Bo = l-6-J2nl (29) 



k 



where n^^ = Np{e^^), Np{e) is the Fermi distribution function, the spectra e^'^k ^^^ deter- 
mined by (|26D with 



■^0 — C "" ^oto 
Ao + /^o = -C - Coto (30) 

and riafi are densities of condensates, index at the parameters stands for their noninter- 
acting mean-field values. It can be checked numerically that the equations (^9[) do not have 
the solutions with positive n^, rib for any 6, (. 

To improve the behavior of solutions of equations (pOf), we consider the corrections to 
above noninteracting mean-field theory owing to fluctuations of Q, V (remember that the 
fluctuations of a were excluded by gauge transformation). Introducing Fourier components 
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Qq,5 = E Q^,^+Se'''''' 

i 



and taking the functional derivatives of action (^51), we obtain for the following expressions 



SJ^tk + sx 



1 nia+nV -^ 



+ 7- 



'^k + ^k+q 



T Z^„ - £^+q + e^ (Zo^n - ^k+q + ^k) ^ , 



Gq (q,^^n)e 



Jk(<5-<5') 



T 



<(! + <) 



(cLcka) = ^k + E 



q,ja;n 



1 72^(1 -n^) 

Tz^„-4 + £tq 



<-q + < 



UC^n 



££ + £ 



k-qJ 



^55' 



G^^(q,za;„)e 



»(q-k)(5-<5') 



(5^tk + '5A 

-<(i-0 



T 



(aj,ak) = < + E 



q,«<^„ 



<(i + 



'^k + ^i^-q 



T iujn -el + ei_^ 



\iuo„ - et + e 



k-qJ 



Gf(q,iu,„)e--<'-''' 



*^*- + ^' + ^„i;(i + „i;) 



T 



(/.I/k) = < + E 

q,ia;„ 



4(1 - 4) 



ni + n" 



k+q 



5Stk + (5A + (5/i 






[lUOr. 



'k+q 



+ £ 



6^2 



G^^(q,^^„)e 



J(k+q)(5-5') 



(31) 



and 



(^1^.) 



(32) 



- E [nM'(q.^^n)-n^y(q,^o;„) 

where (5JF, (5^, SC,SB, SX and (5/i denote the corrections to corresponding quantities of nonin- 
teracting mean-field theory owing to interaction with the fields Q and V; the Green functions 
Gqp are defined by 

Gg'(q,r) = (T[Q,,,(r)Qq,5'(0)]) 
Gf(q,r) = (T[P,,,(r)Pq,,,(0)]) 



and can be found by inverting corresponding matrix in (^). To keep the spectrum of the 

bosons a, b gapless we choose 

1 



6X= E 



q,ia;„ 



lUJr, 



G'^{ci,tuj^^)-6J^to 



SX + 5n= E 



q,iaj„ ^^n ^ — q 



G'j^{ci,tUJn)-dCto 



(33) 



The expressions ( pi]) and (0) together with the self-consistent equations (0), ( ]28|) give the 
parameters of supersymmetric mean-field theory. 

V. CONCLUSION 

In this paper we have obtained the supersymmetric action of t — J model, which inter- 
polates between slave-fermion and slave-boson ones. We have also developed a mean-field 
approximation for this action, which takes into account the fluctuations owing to the in- 
teraction with gauge fields. Physically, these fluctuations correspond to fluctuations of the 
total boson (fermion) number in the system, since only total number of bosons and fermions 
is conserved in the supersymmetric representation [Q. In other words, these interactions 
describe virtual boson-fermion transitions within the supersymmetric fermion-boson repre- 
sentation. The account of these virtual transitions can give in principle a possibility to 
investigate the crossover from the saturated ferromagnetic state into nonmagnetic one and. 
After a generalization to finite U, the crossovers from antiferromagnetic state to the states 
which are described well by slave-boson representation, e.g. pseudo-gap and superconduct- 
ing one, could be described. The application of the results of the present paper to solving 
these problems is the aim of future work. 

VI. APPENDIX. GAUGE TRANSFORMATION OF MEAN-FIELD 
HAMILTONIAN PARAMETERS AND GREEN FUNCTIONS 

In this Appendix we consider the transformation laws for different quantities under the 
gauge transformation (|2^) . We treat rji^r) as a sort of fluctuations (see, e.g. [l^) with zero 
average and the pair correlation function 

/(q, uj) = (r/(q, uj)r]{ci, uj)) (34) 

Calculating the contributions up to quadratic terms in rj, we obtain following transformation 
laws of the Hamiltonian parameters: 

10 



A ^ A - X! (2C - ^^n)/(q, iuJn) 

A^A+ E (^ - •^iq/^0)/(q, lOJn) 
q,ia;„ 

B^B- Y.{B-CtJto)f{q,iu;n) 



C ^ C + ^ (C - i3tq/^o)/(q, ^0On) (35) 



q,2aj„ 



The gauge fields are transformed as 

Qij -^ Qij + Ar]i - Tr]j 

'Pij -^ ^i3 + Cr]i - Br]j (36) 

Using (|36D one can obtain tlie transformation laws for tfie gauge fields Green functions, e.g., 

^^.(q, i^n) -* G'„(q, iun) + (2C - iujnff{(\, iUn) (37) 

Gg'(q, tuJn) ^ GQ(q, «a;„) + (^ - :Fe^'){A - J'e'^'^fiq, tco„) 
Gfici,tco^) ^ GQ(q,«a;„) + (C - Se^'i^)(C - Be'^'')fici,iuj^) 



Tlie transformation laws for other quantities can be obtained by combining the results (^ 
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